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The effect of orthorhombicity on Josephson tunneling in high T c superconductors such as YBCO 
is studied for both single crystals and highly twinned crystals. It is shown that experiments on highly 
twinned crystals experimentally determine the symmetry of the superconducting twin boundaries 
(which can be either even or odd with respect to a reflection in the twinning plane). Conversely, 
Josephson experiments on highly twinned crystals can not experimentally determine whether the 
superconductivity is predominantly s-wave or predominantly d-wave. The direct experimental de- 
termination of the order-parameter symmetry by Josephson tunneling in YBCO thus comes from 
the relatively few experiments which have been carried out on untwinned single crystals. 
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Introduction. Although many of the currently studied 
high-temperature superconductors are orthorhombic (in- 
cluding YBa 2 Cu 3 7 -5 (YBCO) and Bi 2 Sr 2 CaCu 2 8+( 5, 
jlj) they are often regarded as being approximately 
tetragonal. Thus, in the current debate concerning the 
symmetry of the order parameter, the attention has been 
focussed on attempting to distinguish between order pa- 
rameters having s, d x 2_ y 2 or possibly other symmetries 
(see for a recent review). As a result of the orthorhom- 
bicity of the high T c superconductors, however, the s and 
d x 2_ y 2 order parameters are coupled This note dis- 

cusses some of the consequences of orthorhombicity and 
of this coupling in the interpretation of Josephson tun- 
neling experiments, and focuses on what can be directly 
measured experimentally without relying on microscopic 
theories or models of the superconducting state and the 
twin boundaries. 

Experimental studies of the Josephson effect have been 
suggested [|||] as a method of establishing the order pa- 
rameter symmetry in the case of unconventional super- 
conductivity. Measurements of the basal-plane Joseph- 
son currents in YBCO |7|-[TT| have found these currents 
to be characteristic of d-wave symmetry. On the other 
hand, a measurement jl2) of the c-axis Josephson current 
raised questions about the d-wave interpretation. 

Below, I study the Josephson currents from an or- 
thorhombic superconductor such as YBCO to an s-wave 
superconductor such as Pb for both single crystal and 
highly twinned (all twins occurring with equal weight) 
samples. The superconducting state of YBCO can be 
shown to be either even or odd with respect to reflection 
in a twinning plane, and it is demonstrated that Joseph- 
son tunneling experiments on highly twinned samples de- 
termine this symmetry. 

Contrary to current views, however, Josephson exper- 
iments on highly twinned samples can not determine 
whether the superconductivity is predominantly s- or d- 
wave. The Josephson current from a single twin con- 
tains both s- and d-wave components (as it must be- 
cause a single twin is orthorhombic). The averaging over 
all twins, however, effectively imposes tetragonal sym- 
metry and results in a Josephson current which depends 
on only a single symmetry component (either s- or d- 
wave); which component is selected depends on the su- 
perconducting twin-boundary symmetry. Josephson ex- 
periments on highly twinned samples will thus measure 
the magnitude of the selected component, but will give 
no information on the magnitude of the other component. 
In order to determine whether the symmetry of the or- 
der parameter is predominately d-wave or predominately 
s-wave, it is essential to measure the relative magnitude 
of the two components (which can be done on untwinncd 
single crystals, but not on twinned crystals). 

In view of the above, the direct experimental determi- 
nation of the order-parameter symmetry by Josephson 
tunneling in YBCO thus comes from the relatively few 



experiments [pxj|j7|] which have been carried out on un- 
twinned single crystals. 

Mixed s- and d-wave superconductivity. The Ginzberg- 
Landau free energy density describing coupled s-wave 
and d x 2_ y 2 order parameters (tp s and ip d , respectively) 
is 

F = a s |^ s | 2 + a d \4>d\ 2 

+ (-imi^i 2 + /? 2 iv s i 2 k^ + ma 

here e = 1 , 2 refers to one of the two possible twin ori- 
entations which occurs in the YBCO materials. Except 
for the factors of (— 1) £ , which are important below, this 
free energy has been given in [0 . I assume that the twin 
boundary between two twins (as in Fig. 1) is a plane of 
reflection symmetry of the underlying crystal lattice. If 
the d- and s-wave order parameters describing supercon- 
ductivity in twin 1 are (ip d ,ift s ), then a reflection in the 
twinning plane gives a superconducting state of twin 2 
described by order parameters (—ip d ,ip s ) and having the 
same free energy as the corresponding state of twin 1; 
hence the factor (— l) e which appears at various places 
in the free energy. 

By definition, d-wave superconductivity will be mod- 
elled by the assumption that a d — a' d (T — Tdo), while 
a s is taken to be positive. (Conversely, s-wave supercon- 
ductivity will be modelled by taking a s — a' s (T — T cs ) 
with a d > 0.) For d-wave superconductivity, minimizing 
F with respect to tp* in zero external magnetic field gives 
tp s correct to terms linear in ip d as 

i> a = (-l) e+1 (a/a s )ip d 

Also \ip d \ 2 — ct' d (T d — T)/(3 d ; here (3 d is a constant ap- 
proximately equal to (3 d if a is small, and T d differs from 
T d Q when a is non zero. 

Since ip s can be found as an expansion in powers of ip d 
(and this is true for both s- and d-wave superconductivity 
as these have been defined above) , ip s can be substituted 
for in terms of ip d only (the coupling to ip s nevertheless 
being implicitly taken into account). 

The nature of the superconducting twin boundary be- 
tween twins 1 and 2 will be determined by the interface 
free energy per unit area 

Fu = B{ip* dl i) d2 + ipdiiph)' 

where, as just explained, ip s is implicitly taken into ac- 
count. If B < (or B > 0), the free energy is minimized 
by taking ip dl = tp d2 (or i/'di = —ipd2), and the super- 
conducting state will be odd (or even) under reflection 
in the twin boundary, respectively. Note that if the twin 
boundary is odd, ip d i = ip d 2, but that because of the fac- 
tor (— l) e in the above free energy, ip s i = —ip S 2- Thus, 
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calling a twin boundary a normal or a tt boundary ac- 
cording to whether or not the order parameter changes 
sign could be ambiguous, and it is better to identify the 
type of twin boundary by its reflection symmetry. 

Basal-plane Josephson currents. The twin boundaries 
in YBCO are normal to the [110] and [110] directions. 
A YBCO crystal will consist of some regions where or- 
thorhombic twins are separated by [110] twin boundaries, 
and other regions where orthorhombic twins are sepa- 
rated by [110] twin boundaries. For want of a better 
name, I will use the word "macrodomain" to refer to a re- 
gion where the twin boundaries all have the same orienta- 
tion. Thus twinned YBCO can be considered to be made 
up of [110] macrodomains and [ 110] macrodomains. I be- 
gin by discussing a [110] macrodomain, which contains 
[ 110] twin boundaries and two types of twins (twin 1 and 
twin 2) as shown in Fig.l. 

Consider a Josephson junction between twin 1 of 
YBCO and Pb. The normal n (directed from YBCO 
to Pb) to the plane of the junction is taken to lie in the 
YBCO a-b plane and to make an angle with the a-axis 
of twin 1 (see Fig.l). The Josephson current from twin 
1 to Pb is 

3i(0) = ig(9)(tp dl ip Pb - ipdiipU)- 

Now write 

g(0) = d(0) + s(0) 

where 

d(0) = \g(0)-g(6 + */2)]/2, 

s(0) = [g(0) + g(0 + -K/2)]/2. 
The constraints of orthorhombic symmetry require 
g(0)=g(~e), g(0 + n)=g(0), 

d{9 + ir/2) = -d(9), s(0 + ir/2) = s(9). 

It is clear from these symmetry properties of d(9) and 
s(9) (and this will appear in more detail below) that these 
functions will exhibit characteristic d x 2_ !/ 2-wave and s- 
wave behavior, respectively, in Josephson tunneling ex- 
periments; the above separation of g(9) into d(6) and 
s(9) is thus a fundamental step in the interpretation of 
the basel plane Josephson tunneling from an orthorhom- 
bic superconductor to an s-wave superconductor in terms 
of d- and s-wave components. The following Fourier se- 
ries representations of the functions g(9), d(9) and s(9) 
incorporate the above-mentioned symmetries: 

g(9) = S g n cos(2n0), 
d(9) = S d 2+in cos[(2 + 4n)0], 



s(9) = S s 4 „cos(4n6>). 

In all cases the sums are over all integral n from zero to 
infinity. 

The Josephson current from twin 2 to Pb in the direc- 
tion 9 can be obtained from that for twin 1 by reflecting 
in the twin- twin interface (see Fig. 1). Hence 

32(9) = -igfa + S - 9)(r d2 ^Pb ~ Mh)- 

For the purpose of interpreting corner squid experi- 
ments, the current along two orthogonal directions in the 
basal plane (which will be called the x and y directions) 
will be calculated; the x and y directions are chosen to 
be at angles of 6 = 6i + S/2 and 9 = 9 X + (tt + 5)/2 
from the ai direction (see Fig. 1). The angle 9\ (just 
defined in terms of 9 by the preceding equations) is used 
rather than 9 to define the x and y directions because 
9\ = is a symmetrical orientation of x and y relative to 
the two twin orientations. Finally, the expressions given 
below for the Josephson currents will be for currents av- 
eraged with equal weights over the two twins in a [110] 
macrodomain, and will be given in the form 

javg,a = C a (9 1 )\lp d 1p Pb \sin((p d - <j> Pb ) 

where a is a; or y, ip d = \ip d \exp(i<fi d ), and ip Pb = 
\ipP b \exp(i<j>p b ). 

For the case of even twin- twin interfaces, ipdi = 
-ipd2 = tpd, and 

cx(0i) = +[d(9+) - d(9^)] + [s(9+) + s(9^)], 

CyWt) = ~[d(0+) - d(6i)\ + [8(0+) + s(6i)l 

where 9+ = 9\ ±5/2. For odd twin- twin interfaces, for 
which tpdi = 4>d2 = ipd, 

c x (0i) = +[d(9+) + d(0f)] + [s(9+) - 8(0^)], 

cyfa) - -[d(6+) + d(9^)} + [s(6+) s(0i)}. 

The above discussion (which assumed that the entire 
crystal was made up of a single [110] macrodomain) 
can be extended to the case where both [110] and [110] 
macrodomains are present. By assuming that the [110] 
and [110] twin boundaries are orthogonal, the expres- 
sions for c x and c y for the [110] macrodomain can be 
shown to be identical to those for the [ 110] macrodomain, 
except that 8 is replaced by its negative. Thus, if an av- 
erage over equally weighted macrodomains is taken, and 
the phase of ipd is assumed to be the same in both types 
of macrodomains, the terms in c x and c y which are odd 
in (5 drop out, and results characteristic of perfect tetrag- 
onal symmetry are obtained. For even twin boundaries, 
Josephson currents with s-wave symmetry are obtained, 
i.e. 

c*(0i) - Cy(0l) = 8(0! + S/2) + s(-0! + 6/2) 
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where an overbar on c indicates that an average over 
the two macrodomains has been taken (in addition to 
the previous average over the two twin types for each 
macrodomain). For odd twin boundaries, currents with 
d x 2_ y 2-w&ve symmetry are obtained, i.e. 

c*(0i) = -c v (9i) = d{9 l + 5/2) + d(-9i + 8/2). 

On the other hand, if ipd changes sign on going from the 
[110] macrodomain to the [TlO] macrodomain, then even 
twin boundaries yield d xy -wave currents with 

cx(0i) = -cy(9i) = d{Qx + 5/2) - d(-9 x + 5/2). 

while odd twin boundaries yield L z -wave (where L z is the 
z-component of an angular momentum) currents with 

c x {6i) = cy{fix) = 8(0! + 5/2) - s{-6i + 5/2). 

Now suppose that the basal-plane Josephson currents 
in YBCO, when determined on samples in which all twins 
and macrodomains are present with equal weights, are 
characteristic of d x i_ v i -wave symmetry, as is suggested 
by experiment Q PJll| . From above this implies that the 
twin boundaries are odd and tpd does not change sign 
between [110] and [110] macrodomains. 

It is to be emphasized that the experiments on heav- 
ily twinned samples do not determine experimentally 
whether the superconductivity in orthorhombic YBCO 
is predominately s-wave or predominantely d x 2_ J/ 2-wave. 
Above, the function g{9) characterizing the basal plane 
Josephson currents was separated into two components, 
a d-wave component d(9) and an s-wave component s(9). 
However, as shown above, in the case of odd twin bound- 
aries and ipd not changing sign between macrodomains, 
the average contribution of the function s(9) to the 
Josephson currents is zero. Thus, the experiments mea- 
sure only the d-wave component d(9), independently of 
whether d(9) is much greater than or much less than s(9). 

Josephson tunneling experiments on untwinned single 
crystals [|7|JlC|l can, of course, classify the superconductiv- 
ity as being predominantly s- or d x 2_ y 2-wave. 

Finally, note that calling the superconductivity pre- 
dominantly s- or d-wave according to the relative aver- 
age magnitudes of s(9) and d(9), while an appropriate 
definition from the point of view of tunneling experi- 
ments, is only one possible definition (and one in which 
the anisotropy of the tunneling matrix elements will play 
a role) . Another definition which is not necessarily equiv- 
alent could be based on the relative magnitudes of the s- 
and d-wave contributions to the gap function. 

Josephson currents in the c direction. Symmetry con- 
siderations show that the c-axis Josephson currents for 
twins 1 and 2 (e.g. see Fig. 1) have the form 

jze = C Z (-1) £ (V4V , P6 - ^delppb)- 

For odd twin boundaries, ipdi = ipd2, and the Josephson 
current averages to zero when averaging over an equal 



weighting of twins 1 and 2. However, in the c-axis mea- 
surements of Josephson tunneling [Q, the Josephson 
current is found not to average to zero. Thus, these ex- 
periments would seem to imply that the twin boundaries 
are even, in contradiction to the result implied by the 
above interpretation of the basal-plane Josephson cur- 
rent experiments. 

Conclusions. Josephson tunneling experiments mea- 
suring Josephson currents from an orthorhombic super- 
conductor such as YBCO into an s-wave superconductor 
such as Pb, and carried out on samples in which all twins 
and macrodomains have equal weights, experimentally 
determine the twin-boundary reflection symmetry of the 
superconducting state. In the case of YBCO, the basal- 
plane experiments [Q,^| characteristic of d-wave symme- 
try suggest that the superconducting state has odd re- 
flection symmetry with respect to the twinning plane, 
whereas a c-axis Josephson experiment jijj suggests that 
the superconducting state has even twinning-plane sym- 
metry. 

A similar analysis carried out for YBCO to YBCO 
Josephson tunneling shows that here also the experi- 
ments on twinned crystals experimentally determine the 
twin boundary symmetry; the results of [QJ|,|l^] suggest 
that the twin boundaries have even symmetry, in con- 
trast to the results of Jl3| ] which suggest even-symmetry 
twin boundaries. 

The above interpretation assumes that, in twinned 
crystals, all twins and macrodomains are present with 
equal weights. If this is not true, the interpretation is 
not valid. For example, if, in the experiment of Sun et 
al. pS^], twins of type 1 occupied an area of the Joseph- 
son junction which was significantly larger than that oc- 
cupied by twins of type 2, the assumption of odd sym- 
metry twin boundaries would be acceptable since there 
would then be only a partial cancellation of the Joseph- 
son currents, and a nonzero Josephson current would be 
observed. The results of the c-axis Josephson experi- 
ments and the basal-plane Josephson experiments could 
then be reconciled. Since this question of the weighting 
of the different twins and macrodomains is of crucial im- 
portance to the interpretation, it would be of interest to 
try to determine this weighting experimentally, for exam- 
ple by microstructural studies of the Josephson junctions 
used for both the c-axis and the basal-plane experiments. 

Another important conclusion is that, contrary to cur- 
rently accepted ideas, Josephson tunneling experiments 
on orthorhombic YBCO crystals containing all twins and 
macrodomains with equal weights can not directly ex- 
perimentally determine whether the superconductivity is 
predominately d- or s-wave (but can, as described above, 
experimentally determine the nature of the twin bound- 
aries). Our direct experimental knowledge of the order 
parameter symmetry as determined by Josephson exper- 
iments thus comes from experiments performed on un- 
twinned single crystals 
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Fig. 1 (a) Shows a [110] twin boundary(DB), the two 
twins with the directions of their a and b axes, and 
the interface between YBCO and Pb (ABC) with 
its normal n. (b)Definition of the angles between 
ai, a 2 , and n. 
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